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We use Bloch wavefunctions in a new way to study a fluid interfacial problem and find both linear 
oscillatory and nonoscillatory instabilities of the surface deformation. Underlying periodic flows, 
such as those arising in Rayleigh-Benard cells, Langmuir circulation, and solar convection rolls are 
treated as regions of varying surface shear which scatter and refract surface capillary-gravity waves. 
We find that Bloch wavefunction decompositions of the surface deformation n(r) and velocity po- 
tential tp(r, z) results in a nonhermitian secular matrix with an associated band structure that gives 
rise to extremely rich surface instabilities (complex frequencies) . These are predominantly enhanced 
at certain Brillouin zone edges, particularly near Bragg planes corresponding to the periodicity de- 
termined by converging or diverging surface flows. The instabilities persist even when the dynamical 
effects of the upper fluid are neglected, in contrast to the uniform shear Kelvin-Helmholtz instability. 
The periodic flows can also couple with uniform shear and suppress standard Kelvin-Helmholtz (K-H) 
instabilities. 
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Introduction - Regular fluid flow patterns are ubiqui- 
tous . Periodic flows are typically the result of thermal 
and/or dynamical instabilities and can exist in parameter 
regimes prior to chaotic or turbulent flows. An example 
on laboratory scales is the well-known Rayleigh-Benard 
convective instability ||]. Heating of a fluid layer from 
below provides buoyancy which induces the instability. 
As the instability develops, regular convective patterns 
in the form of rectangles and hexagons appear |^|| . 

Another physical example of nearly periodic flows 
with a free surface occurs in Langmuir circulation (LC) 
windrows where wind stresses, turbulent stresses, and 
coriolis forces conspire to form convection rolls in the 
upper ocean 

H|). These nows 

can be spatially periodic 
as often observed when the sea surface is contaminated 
by oil slicks or algae. The LC process is believed to be 
an important mechanism in thermal and chemical mixing 
in the upper ocean. Mechanisms of LC generation may 
involve the interaction of turbulent surface currents with 
surface waves which in turn mediate the wind stresses 
H . Refraction of largely irrotational surface waves from 
the underlying rotational Langmuir circulation rolls are 
thought to contribute to wave breaking and enhanced lo- 
cal wind stresses On an even larger scale are the 
approximately periodic solar convection cells in the so- 
lar convective zone. In this system, magnetic fields will 
likely also affect the surface wave dynamics f|. 

Clearly, surface wave propagation and interfacial sta- 
bility is an important aspect of stratified fluids with wide 
applicability. Oscillatory instabilities of the two-layer 
Benard problem have been studied by considering all cou- 
pled modes, including thermal and Marangoni In the 
limiting one layer Rayleigh-Benard problem, Benguria 
and Depassier even find interfacial oscillatory insta- 
bilities for parameters prior to the onset of the periodic 
roll states. This occurs in the fixed lower temperature 



and fixed upper thermal flux ensemble. In this Letter, 
we assume a pre-existing periodic flow and explore its 
dynamical effects on interfacial instabilities. By consid- 
ering surface wave motions as a perturbation to the un- 
derlying mean flows, we treat the surface waves as being 
reflected or refracted much like wave scattering in opti- 
cal, acoustic, or electronic solid state physics, where band 
structures have been calculated using linear eigenvalue 
analysis in various systems A related problem of pe- 
riodic surface wave scatterers such as thin ice floes has 
been treated with similar methods || . Here the underly- 
ing flow is periodic so we use Bloch functions to describe 
the surface displacements and dynamic fluid pressure and 
derive a generalized quadratic eigenvalue equation with a 
nonhermitian operator. In addition to nontrivial "band 
structure" we find a complex eigenvalue spectrum corre- 
sponding to linear instabilities with rich behavior. We 
show the existence of interfacial instabilities in the one 
layer problem, as well as modifications to the stability 
regions of uniform shear in the two-layer problem (K-H 
instability) when an additional periodic flow is present. 




FIG. 1. Schematic of problem geometry. Two ideal fluids 
are separated by an interface with surface tension a under the 
action of gravity. The densities, velocities, and depths of the 
upper and lower fluids are , U ± , and respectively. 
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Formulation - Consider two ideal fluids separated by 
a flat wave-undisturbed interface at z = as depicted in 
Fig. 1. Although the underlying flow field is typically 
rotational, we assume the flow associated with the im- 
posed surface waves (externally excited) are irrotational. 
The displacements caused by the surface waves are as- 
sumed small and do not appreciably affect the under- 
lying flow. This requirement, and linearization, implies 
?y <C A, a where 77, A, and a are typical surface amplitudes, 
wavelengths, and flow field periodicities respectively. 

The fluid velocities above and below the interface are 
= U + Vtp^ where U is the periodic flow field (gen- 
erated for example, by the Rayleigh-Benard instabilities 
0, Langmuir circulation mechanisms or electrohy- 
drodynamic effects jlO|]) which satisfies V-U = 0, and 
ip is the velocity potential for the irrotational capillary- 
gravity waves, and z = i/i^ is the position of the impen- 
etrable top and bottom. For concreteness, consider the 
Rayleigh-Benard system where instabilities to periodic 
flows in the Boussinesq approximation with a free surface 
arises when Ra = agh 3 AT/(vK) > Ra* , (Ra* ~ 1100 for 
a free, tensionless surface) where a, g, AT, v, and n are 
the thermal expansion coefficient, gravitational accelera- 
tion, temperature difference T(z = 0) — T(z = ~h), kine- 
matic viscosity, and thermal conductivity, respectively. 
The ideal fluid approximation is valid only for surface 
waves which are not significantly damped over many pe- 
riods of the underlying flow, i.e., for capillary waves, the 
attenuation length k^ 1 = 3a / '{Av p~ lu) ^> a. We also as- 
sume in this system a large Prandtl number such that 
the Reynolds number, Re ~ yTla/Pr is not too large as 
to render the steady periodic flows unstable. 

Incompressibility demands V • v = Aip = (Aj_ + 
d 2 )^ = 0, where A = A j_ + d 2 is the three dimensional 
Laplacian. The boundaries at z = ±/i ± are assumed flat 
and impenetrable. Therefore, the linearized kinematic 
boundary conditions at z = ±/i and the interface at 
z ~ are 9 2 <p ± (af, z = =ph) = and 



harmonic source can be found by superposing the solu- 
tions of many frequency components. Consider general 
solutions to (ys* and r](x), 



± / -* \ ± iq r C0S ^ 1 l(^ ± T z ) 

ip (r, z ) — > iprze q : — - 

y V ' ; 4^ 9 coshg/i± 



and 



77(F) = 



(3) 



(4) 



where q = q± lies in the surface plane, and q = 
lim e ^o y/ q 2 + e 2 . Equation (|J) automatically satisfies 
\7 2 ip — as well as the impenetrable bottom boundary 
condition. Therefore by using Eqn. (^), the problem is 
reduced to that of simultaneously solving equations (|l|) 
and (||) with the unknowns tp^ and r\q. A Fourier decom- 
position for the periodic flows is 



z )e 



G-r 



(5) 



where G is the reciprocal lattice vectors appropriate for 
the underlying flow periodicity. With the form (|^), 
the velocity potential <p will furthermore be restricted 
to forms satisfying Bloch's theorem tp(f + a,z) = 
e lq ' r f(f, z) where f(r) is a function periodic with respect 
to translations of a. 

We wish to find the complex dispersion relation be- 
tween lu and wavevector q. Substituting Eqns. (|^) and 
(El) into (TO) and (0), we obtain 
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and 



± 



qtanhqh^ 



(G) 



d t r l (i f ) + U ± (r)-V ± r)(r) 



lim d z tp ± {f, z) (1) 



respectively, where fj_ = r = (x,y). The linearized dy- 
namical boundary condition at z ~ is found by balanc- 
ing z— component stresses from the dynamical pressure 
with that from gravity and surface tension; 



lim \d t p ± + U ± d l >p ± ] = A 



-LV - grj, 



(2) 



where we have for simplicity assumed a constant surface 
tension and neglected any possible Marangoni effects that 
may arise when surfactants are convected along the sur- 
face by U(f, z = 0). Spatially varying surface properties 
such as tension or bending rigidity can be treated without 
difficulty [p|Jll|. Henceforth, we assume that all quanti- 
ties vary as e~ tuJt . Wave evolution due to a non-time 



icu<p~ — % 



Y^U-(G)-(q-G)p ? _ d 



T+9) Vq 



(7) 



respectively. Shifting the the wavevector into the first 
Brillouin zone, and substituting (|^) into (0), we find a 
generalized eigenvalue problem 



^2 (Ac^ 2 + Buj + C) rj^(G') = 



(8) 



where the matrices in the limit of unbounded upper fluid 
(h + = 00, h = h~) are given by 
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A(G, G') = {s tanh \q - G\h + 1)S S s , 



(9) 



B(G, G') = -2sU+-{q-G) t&ah\q-G\h5 3s 
\q-G\i&nh\q-G\h ' 



1 



|g-G'|tanh|g-G'|/i 



U-{G-G')-{q-G') 



(10) 



G(G, G') = (s|t/ +-(g- G)| 2 tanh |g - G> - f>|(G)) 
+ ^ |g-G*|tanri|g-G|/i 



-U-(G-G")-tf-G") 
q-G"\tanh\q-G"\h 

x IJ-{G" -G')-{q-G') 



where s = p + / p , and 



(11) 



fi|(G) = y—W- G\ 3 + (1 - s) 9 |<f- G|) tanh \q - G\h 

(12) 

We arrive at an eigenvalue problem that can be solved 
by standard means p2| , 



1 

-A - l C A _1 B 



Vjr 



(13) 



where tp = ujff. However, the above matrix is gener- 
ally nonhermitian and the corresponding spectrum lu is 
complex. The standard criteria for the Kelvin- Hclmholtz 
propagating wave instability is recovered when U ± = U Q 
are uniform jl3| : 

(uj-U^-kf + steinhkh(uj-U+-k) 2 -ftfc(0) = 0. (14) 

Note that when dynamical effects of the upper fluid are 
neglected, (s = 0), the roots of Eqn. (14) are real, and 
no linear instability exists. 

Results and Discussion - Distances, frequencies w, and 
velocities U will be normalized and measured in units 
of a, yj g/a, and ^fga respectively. First consider only 
periodic flow in the lower fluid with no uniform com- 
ponent, J7 ± (G = 0) = 0. Only the velocity at the sur- 
face U~ (r, z — 0) will influence surface wave propagation. 
For simplicity we only analyze one-dimensional rolls de- 
scribed by the approximate function, 



U (r,z = 0) = U (z)xcos 



2n 



(15) 



consequently, U (G) = U (0)x/2tt for G = ±2nx/a, 
and zero otherwise. This choice of phase for U~ (r , z = 0) 



implies one converging and one diverging surface flow 
region per unit cell and also simplifies the computation by 
making all elements of the matrix in Eq. ( [Df ) real. Also, 
we choose small surface tensions so that the periodic flow 
surface velocities will be nonnegligible compared to water 
wave group velocities. 
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FIG. 2. Triptych depicting the band structure for periodic 
flow in the lower fluid (s = 0, h = 2.0, Bo" 1 = 0.01). The 
central panel shows the dispersion relation in the q x direction 
in the reduced zone scheme, (a). U~(2n/a) = 0.20; (b). 
U~(2ir/a) = 0.75. 

We find the spectrum of (jl^) using standard meth- 
ods jl2| and thus obtain the band structure of capillary- 
gravity waves on surfaces with underlying periodic flows 
U~ . The central panels in Figures |^ show the real (solid 
lines) and imaginary (height of hatched regions centered 
about Re{w}) parts of ui(q x ,0). We use the notation q (k) 
to denote quantities plotted in the reduced(extended) 
zone scheme. The side panels show w(0, k y ) for = 0, 
h = 2.0, s = 0, and the inverse Bond number Bo" 1 = 
a/(p-a 2 g) = 0.01. U-(±2n/a) = 0.2, 0.75 are shown in 
(a), and (b). respectively. The "band structures" shown 
in 2(a). contain branch cuts at certain q x satisfying the 
Bragg scattering condition. There are open band gaps 
at q x = 0, which decrease at larger u, similar to elec- 
tronic and acoustic wave propagation in periodic media; 
the gaps normally found at q x — ±ir are collapsed due 
to the converging flow in each unit cell and are degener- 
ate down to a smaller value of q x . Under this periodic 
flow, growing surface modes arise near q x ~ ±7T when 
real degenerate roots of Eqn. (jl^) split into complex 
conjugate pairs, even though for the one fluid problem, 
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uniform stream flows are stable according to ( |l4| ) . In Fig 
§(b), \U~ (±2ir / a)\ — 0.75 where the first gap at q x = 
has also merged. Note that unstable modes are associ- 
ated with both standing and low group velocity travelling 
waves since they appear predominantly near q x = 0, it. 
These may be damped or saturated by viscous dissipa- 
tion, Im{w} ~ 2v\q- G\ 2 . For tf-(2ir/a) < 0.3641, or 
higher values of surface tension, the band structure quali- 
tatively resembles that of Fig. 2(a). At intermediate val- 
ues of U~(2n/a) (not shown), we find unstable nonoscil- 
latory modes; i.e. Re{tu(q x , k y )} = 0, Im{oj(q x ,k y )} < 
0. When 0.417 > tf-(2n/a) > 0.3641, the lowest 
two branches collapse such that unstable zero frequency 
modes proliferate and fill the whole zone. These growing 
modes will be temporally saturated by viscosity and/or 
nonlinear effects, and result in static surface deforma- 
tions similar to the Reynolds ridge juj. Upon further 
increasing U^(2n/a) > 0.417, the zero frequency grow- 
ing modes develop a finite frequency and stabilize near 
q x ~ 0. When U~{2n/a) ~ 0.428, the first gap at q x = 
between the stable modes collapse. The qualitative be- 
havior described above continually repeats upon increas- 
ing U~ (2tt I a). Note that for higher U~(2-k/o) there are 
specific directions where oscillatory and nonoscillatory 
(Re{u)(q x , k y )} — > 0) instabilities also arise. 




FIG. 3. Dispersion relation with and without an under- 
lying periodic flow with a uniform (Uq = 2.0) overlying 
shear. Here, h = oo, Bo -1 = 0.03, and s = 0.1. As a 
visual aid, we have imposed a k x dependent Galilean shift 
A(7o = (Uq + UgS tanh k x h)/(l + st&nhk x h). The pure 
shear dispersion {U~ = 0) is shown by the dotted lines, 
with the imaginary parts of the unstable growing modes de- 
picted by the heights of the hatched regions. The solid line 
(u) < 0) defines the dispersion relation when a periodic flow 
(U~(G 7^ 0) = 0.2) is imposed. Note the changes to the 
instability regimes. 

The band structures shown in Figures 2 are only quan- 
titatively altered (tilted) when a small uniform shear 
is imposed in addition to the periodic flow. However, 
for uniform shear (s = 0.1, h = oo,Uq — 2.6, Bo -1 
= 0.1) flows with an existing K-H instability (with un- 



stable wavevectors near k x ~ tt) as determined by Eq. 
(|t|) , the effect of an additional periodic flow in the lower 
fluid enhances these instabilities. However, when the 
uniform shear instabilities span k x ~ 27T (s = 0.1, h = 
oo, Uq = 2.0, Bo -1 = 0.03) where there is an open gap, 
the K-H instability can be suppressed by a periodic flow 
U~(2Tr/a) = 0.20, although more instabilities arise at 
higher k x as well. These effects are shown using the ex- 
tended zone scheme in Figure 3. The structure of the 
dispersion relation is rather sensitive to the amount of 
underlying periodic flow and can change drastically with 
variation in any of the parameters. 

Qualitatively, the upward convex tongue of the K- 
H instability in the k x — plane is modified when 
U~(G ^ 0) is added. For instabilities straddling the 
open gap Bragg planes, this tongue is shifted, with parts 
shifted upwards (resulting in destabilization) and down- 
wards (resulting in stabilization). The discontinuity in 
growth rate is clearly shown in Fig. 3 near k x ~ 27r. 
Moreover, other smaller tongues develop at higher q x 
(k x ~ 37T in Fig. 3) corresponding to the closed gap in- 
stabilities that would occur in exclusively periodic flows 
(Fig. 2). The growth rates within the instability regions 
change as well as the delimiting instability regions. Qual- 
itatively, these effects on K-H instabilities result from 
suppression of the travelling unstable waves in K-H shear 
flow, due to coupling to the standing waves of periodic 
flow, arising in the last term in C in Eq. (|i"l|). 

We have shown that periodic flows couple different re- 
ciprocal wavevectors together and affect interfacial sta- 
bility in a nontrivial manner. In particular, by solving 
the matrix equations, we find, remarkably, that an in- 
terface overlying periodic flows is generally linearly un- 
stable, even if p + = 0. The linear instabilities described 
would be saturated by both viscous and nonlinear effects. 
Regions of instability corresponding to finite frequency 
modes with vanishing group velocities, dissipation due 
to viscosity will result in a standing oscillating mode. 

The analogies between wave propagation in periodic 
media and surface wave propagation in periodic shear 
flow configurations can be extended to consider more 
complicated periodic flow structures such as rectangu- 
lar and hexagonal patterns, where rich behavior should 
be expected. This may lead to insights with wide appli- 
cability, from Rayleigh-Benard convection, Langmuir cir- 
culation, solar convection cells, and MHD surface Alfven 
waves in the presence of periodic magnetic fields (in ad- 
dition to flows) || . Furthermore, the influence of defects 
and disorder in the periodic surface flows fl5|| can be con- 
sidered to study surface wave localization in the presence 
of random U(f, 0) (and hence random B, C). For exam- 
ple, methods used to determine the complex spectrum 
density of states of random matrix operators Jl6| should 
yield interesting behavior regarding the sea surface wave 
spectra in the presence of random underlying turbulent 
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surface flows. 
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